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Abstract 

We introduce real second-order freeness in second-order noncommu- 
tative probability spaces. We demonstrate that under this definition, 
three real models of random matrices, namely real Ginibre matrices, 
Gaussian orthogonal matrices, and real Wishart matrices, arc asymp- 
totically second-order free. These ensembles do not satisfy the complex 
definition of second-order freeness satisfied by their complex analogues. 
We use a combinatorial approach to the matrix calculations similar to 
the genus expansion for complex random matrices, but in which nonori- 
entable surfaces appear, demonstrating the commonality between the 
real models and the distinction from their complex analogues, motivat- 
ing this distinct definition. In the real case we find, in addition to the 
terms appearing in the complex case corresponding to annular spoke 
diagrams, an extra set of terms corresponding to annular spoke dia- 
grams in which the two circles of the annulus are oppositely oriented, 
and in which the matrix transpose appears. 



1 Introduction 

In a noncommutative probability space, freeness is an analogue of inde- 
pendence in a classical probability space. Many important random matrix 
models, both real and complex, are asymptotically free; that is, as the size of 
the matrix becomes large, independent matrices satisfy freeness conditions 
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on the expected values of the traces of their products ^23j. No modification 
of the definition is required for real random matrices. 

In addition to considering the asymptotic moments of random matri- 
ces, that is, the expected values of traces of products in the large matrix 
limit, it is also possible to consider the distributions or fluctuations around 
these expected values (central limit-type theorems as opposed to law of 
large numbers-type theorems). Many important matrix models (including 
all of those considered in this paper) are asymptotically Gaussian, so the 
most important quantity of these fluctuations is the appropriately rescaled 
asymptotic covariance of traces. Second-order probability spaces, that is, 
noncommutative probability spaces equipped with a bilinear function mod- 
elling this covariance, were introduced in and further studied in [6|[T7j. 
Also given are definitions for second-order freeness and asymptotic second- 
order freeness, and it is shown that several important complex matrix models 
satisfy this definition. Second-order freeness then has the role for fluctua- 
tions that first-order freeness has for moments. In particular, if random 
matrices A and B are asymptotically second-order free, the asymptotic fluc- 
tuations of ^ + and AB can be calculated from the asymptotic moments 
and fluctuations of A and B. 

Asymptotic second-order freeness as defined in [Tsj is not generally sat- 
isfied by real ensembles of random matrices. If random matrices A^^n and 
Bi^N, k,l = 1, . . . ,p are elements of the algebra generated by a model studied 
in this paper, then the relation satisfied instead is 

lirn^ cov (ti (^ii,Ar • • • Ap^N^ , Tr (^Bi^n ■ ■ ■ -Bp.Af) ) 
p—i p p 
= E n ^ {kNB,..,N) ) + E n ^15^ (tr (i,;vi3j+,^) ) 

fc=Oi=l k=Oi=0 

(where the circle above the random matrix terms indicates that they have 
been centred: X := X — K (ti {X)), cyclically adjacent terms in each entry 
come from algebras generated by independent ensembles, and all indices are 
taken modulo p), while in the limit covariances of cyclically alternating terms 
with different numbers of terms vanish, as in the complex case. (Below we 
will generally suppress the index N.) The first sum on the right-hand side 
appears in the definition of complex second-order freeness, and its summands 
correspond to "spoke diagrams" on an annulus, as in the first row of Figure [T] 
and described in |l8j. The second sum does not appear in the complex 
case. Its summands correspond to the spoke diagrams in the second row of 
Figure [l] in which the two circles of the annulus are oppositely oriented. 
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Figure 1: Spoke diagrams on cycles with three terms. 



We thus propose a definition of real second-order freeness of subalgebras 
of an abstract second-order probability space {A,ipi,{p2) equipped with a 
linear involution reversing the order of multiplication: 

p— 1 p p^i p 

k=Oi=l fc=Oi=l 

(where the aj and bj are centred and from cyclically alternating subalgebras, 
and again all indices are taken modulo p), and similar quantities vanish when 
the number of terms in each argument of ip2 are different. 

In Section [2| we present the notation and definitions we will be using, 
including the relevant existing definitions in free probability, asymptotic 
freeness and second-order probability spaces. In Section [3j we define the 
three matrix models which we will be considering. We find that each model 
can be put into a common form resembling a genus expansion, but in which 
terms corresponding to both orientable and nonorientable surfaces appear. 
We will use this form in the remainder of the paper, in which we will derive 
results which are applicable to each of the three models or any other model 
satisfying the same conditions. In Section |4] we present several combinato- 
rial expressions which are exact for all matrix dimensions N. We derive an 
expression for expected values of products of traces of several independent 
matrices and for the cumulants of such traces, and we find a combinatorial 
characterization of the terms which contribute when the terms are centred. 
In Section [5j we consider asymptotic behaviour. We characterize combi- 
natorially the highest order terms and show that all of the matrix models 
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have second-order limit distributions, meaning that not only the moments 
but also their fluctuations are well-defined asymptotically. In Section [6j 
we show that the matrix ensembles are asymptotically free under the usual 
definition. Using the machinery of [l6], we show that highest order terms 
correspond to a certain type of annular-noncrossing diagram (which we gen- 
eralize to the real-matrix context, in which a larger class of diagrams are 
possible). We show that the highest order terms then correspond to a larger 
set of spoke diagrams, (those for both relative orientations of the two cir- 
cles), and we define asymptotic real second-order freeness as satisfying such 
a relation in Section [7] where we will collect definitions for the real case. 

2 Notation and definitions 
2.1 Combinatorics 

We denote the interval of integers { 1 , . . . , n} by [n], and for m < n the 
interval {m, . . . ,n} by [m, n]. 

If I is a set of integers, then we define — / := {—k : k G /}. We define 
±1 := /U (-/). 

Definition 2.1. For a set /, a partition of I is a set {Vi, . . . , T4} of subsets 
of I called blocks which are disjoint, nonempty, and have union /. We denote 
the set of partitions of / by P (/) , and the set of partitions of [n] by V (n) . 

We define a partial order ^ on "P (n) by letting vr ^ p if each block of vr 
is contained in a block of p. 

For a set I, we denote the largest partition 1/ := {/}. If / = [n], then 
we let In '■= {[n]}- 

We denote the join of two partitions ir, p £ V (n) by vr V /?, where the 
join TT V p £ V (n) is the smallest partition such that vr, p ^ vr V p. 

A partition is a pairing if each of its blocks contains exactly two elements. 
We will denote the set of pairings on a set / by V2 {!) , and the set of pairings 
on [n] by 7^2 (n). If n is odd, then V2 {n) is the empty set, and any sum over 
7^2 (n) is zero. 

We will denote the group of permutations on a finite set / by 5 (/), and 
the group of permutations on [n] by 5„. 

We will often use cycle notation for permutations, in which the permuta- 
tion maps each element to the next element in the cycle (and the last element 
in a cycle to the first). When we multiply permutations, the rightmost acts 
first, then the next to the left, and so on. If we conjugate a permutation vr 
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by another permutation p, then in cycle notation, the resulting permutation 
pirp'^ has the same cycle structure as vr, but with each k replaced by p (k). 

We will also call a permutation a pairing if each of its cycles contains 
exactly two elements. We note that a partition in 1^2 (I) uniquely defines a 
pairing in S (/) and vice versa. 

We denote the number of cycles (orbits) of a permutation tt by # (vr). We 
will use the same notation to denote the number of orbits of any subgroup of 
a permutation group. The subgroup generated by permutations vri, . . . , tt^ 
will be denoted (vri, . . . , vTfc). We can see that # (vr) = # (^~^) and # (it) = 
#(pvrp-i). 

A permutation ir G S (I) is considered to act trivially on a k ^ I. The 
number of orbits of a permutation depends on the implied domain of the 
permutation, which we will state explicitly if it is not clear from the context. 
If vr G 5 (/) and p £ S (J) for / and J disjoint, then # (vrp) = # (vr) + # {p). 

We will often consider the orbits of a permutation as a partition, and 
will use the permutation to denote this partition. 

If J, C I, we say that a permutation or subgroup of S (!) connects 
J and K if it has an orbit which contains at least one element of both J 
and K. We will say that a permutation or subgroup with orbits given by 
partition vr connects the blocks of partition p if vr V p = 1/. 

Definition 2.2. If J C / are sets and vr G S (/), we define the permutation 
induced on J by it, denoted vr| j G S (J), by letting vr| j (k) = vr™ (k), where 
m is the smallest positive integer such that vr™ (k) G J. (In cycle notation, 
this amounts to deleting all elements not in J.) 

While it is not in general true that (vr| j) {p\j) = vrp| j, it is true if at least 
one of vr or p does not connect J and I\J. We can see that (vr|j)~^ = 7r~"'^| j 
and if K C J, vr| = vr|j^. 

Although an induced permutation vr|j is defined only on J, we will find 



it useful in Lemma 4.6 to define vr|j for all A; G / in the same manner: 



vr| J (k) = vr™ (k) where m is the smallest positive integer such that vr™ (k) G 
J even if ^ J. 

Throughout, we will let 5 : k —k. 

Definition 2.3. For a set of nonzero integers / we call a permutation vr G 
5 (±7) a premap if vr (k) = —7r~^ (~^) and no cycle contains both k and —k 
for all k G ±1. We denote the set of premaps on ±1 by PM (±1). 

We note that the first condition is equivalent to dird = vr~^; and given 
the first condition, to show the second, it is sufficient to show that there is 
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no k such that vr (k) = —k (since given a k with vr^™ (k) = —k, we must 
then have vr™" {k) = 7r~™ {—k) = — vr™" (A;), which is impossible, and given 
a k with 7r2"^+i (yfc) = -k, we must have vrvr'" {k) = tt""^ = -vr"^ (fc), 
which is disaUowed) . 

Remark 2.4. The permutation in the above definition is only one component 
of a premap as defined elsewhere, such as in |22j , in which such permuta- 
tions are used to describe surfaced hypergraphs. (Here, a hypergraph is a 
graph with hyperedges, that is, edges which may connect any positive in- 
teger number of vertices, rather than just two, and can be interpreted as 
vertices which must alternate with the other vertices.) In a surfaced hyper- 
graph, not only vertices and faces but also hyperedges must have a cyclic 
order on the graph elements they are connected to. A premap is the appro- 
priate object to describe the faces, hyperedges, or vertices of an unoriented 
surface: the conditions ensure the consistency of the description on an ori- 
entable two-sheeted covering space (see, for example , [l3], pages 234-235, 
for the construction of this two-sheeted cover) . See (7|[l4l[l6l[20|[22] for more 



information on how surfaced hypergraphs can be represented as sets of per- 
mutations. This interpretation is often useful for understanding subsequent 



calculations, but is not necessary to the proofs. See Remark 3.5 



Definition 2.5. We call a cycle of a permutation tt G 5'(ib/) particular if 
its element with smallest absolute value is positive, after jHj. We denote 
by 7r/2 the set of all particular cycles of a premap vr. We will use 7r/2 to 
denote both the set of all elements appearing in particular cycles of tt and 
the permutation on that set by those cycles. 

If vr € S" (/) where / does not contain both k and —k for any k, we denote 
by 7r-(_ the permutation tt considered as a permutation on ±7 (where it acts 
trivially on any element not in J). We define 7r_ := Stij^S. Then vr+vrZ"'^ is a 
premap with # (vr+vrZ^) = 2# (vr). 

We note that the inverse of a premap is a premap. We also note that 
the permutation "fZ^'^"f+ is a premap: 

Lemma 2.6. Let 'j £ S (!) for some set I that does not contain both k and 
-k for any k, and let tt G PM (±7). Then e PM (±7). 

Proof. Firstly, 6^Z^T^"i+5 = 5^Z^55tt55^j^5 = 7^^7r~^7_. 
Secondly, if, for any k, 71^717+ (A;) = —A;, then 

(7+ (^)) = 7- {-k) = -7+ {k) ■ 
But since vr is a premap, this cannot be the case. □ 
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2.2 Random variables and matrices 

Definition 2.7. We define the classical cumulants ki,k2,- - ■ such that the 
nth cumulant A;„ is an n-Iinear function on random variables satisfying the 
moment- cumulant formula, by which they are uniquefy defined: 

E(Xi---X„)= J2 n km{Xi„...,XiJ. (1) 

TreV{n) {ii,...,im}e7r 

In particular, the second cumulant is the covariance: 

k2 {X, Y) = E (XY) - E (X) E {¥) . 

We also note that if any entry Xi of a cumulant kn {Xi, . . . , X„) with n> 1 
is a constant random variable, then the cumulant vanishes. (Briefly, we have 

E{Xi---Xn) = ki{Xi)E{Xi---Xi-iXi+i---Xn) 

= kliXi) Yl n km{Xi„...,XiJ; 

iTer{[n]\{i}) {ii,...,im}eir 

that is, the moment is the sum over terms in which i is in its own block. 
Inductively, we assume that any term in the moment-cumulant formula in 
which Xi is contained in a km with 1 < m < n vanishes, so kn {Xi, . . . , X„) 
also vanishes.) 

We denote the usual trace by Tt (X) := ^^iXa, and the normalized 
trace by tr := ;^Tr, where N is the size of the matrix. 

When we wish to move the subscript of a matrix to its superscript 
(such as when we have indices in the subscript), we will put it in brackets, 
and denote Xk by X^'^\ We will often find it convenient to denote a matrix X 
by X^^^ and its transpose by X^^^\ If we want to combine these notations, 
we will denote X^ by X^''^ and by X^-''\ 

Let / = {ci, . . . , Cm} C =b [n], let mi + ■ • • + = m, and let tt = 
(ci, . . . , CmJ • • • {cmi+-+mk-i+i^ ■■■,Cm) ^ S {I). We denote the trace along 
the cycles of tt: 

Tr^ {Xi, . . .,Xn) 

:= Tr ^X^'^i) • • • x('^i)^ • • • Tr • • • X^'^"'^^ . 

If the domain of vr is restricted to ibJ, where J is a finite set of positive 
integers J = {ji, ■ ■ ■ , jm} with ji < ... < jm, then we may list in the 
arguments of the trace only the elements of J (in order): Tr^ {Xj^^, . . . , Xj^). 
In each case, the domain will be specified if it is not clear from context. The 
normalized trace along a permutation is defined analogously. 
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2.3 Noncommutative probability spaces and freeness 

Definition 2.8. A noncommutative probability space {A,ipi) is an unital 
algebra A equipped with a tracial functional ipi (called an expectation) 
such that (fi {1a) = 1- We will call elements of A (noncommutative) random 
variables. 

When X is a random matrix, we will let tpi {X) = E (tr {X)). 

Definition 2.9. We will say that a random variable a G A is centred if 

(pi (a) = 0. We will denote the centred random variable a, := a — ipi (a). 
In particular, if X is a random matrix, we will denote the centred random 
matrix hj X := X -K (tr (X)). 

Definition 2.10. We say that a product of terms ai,...,ap taken from 
subalgebras Ai, . . . , An C A are alternating if G Ak^ and fci / A;2 / • • • 7^ 
kp, and they are cyclically alternating if, in addition, fcp / fci. Similarly, we 
say that a word w : [p] ^ [C] is alternating if it; (1) ^ w (2) ^ ■ ■ ■ w (p), 
and cyclically alternating if in addition w (p) ^ w (1). 

Definition 2.11. Let Ai, . . . ,An C yl be subalgebras of noncommutative 
probability space A. We say that Ai, . . . , An are free if 

(pi (ai, ...,ap) = 

whenever the are centred and alternating. 

Let (il,S,P) be a probability space. 

Definition 2.12. For each colour c G I for some index set / and for N = 

1,2,.. ., let \x^^^ : n MtvxTV (F)| be a family of matrices. (Here F 

I- J AeAc 

may be either M or C.) We say that the families are asymptotically free if, 

for any integer p > 0, alternating word w : \p\ I, and A^ in the algebra 
generated by the X^^^^ for I < i < p, we have 

JimE(tr(i,...i,))=0. 

Definition 2.13. A second-order probability space {A,ipi,ip2) is a noncom- 
mutative probability space (A, (pi) equipped with a bilinear functional p2 
which is tracial in each argument and such that </?2 {^A, a) = <^2 (a, 1a) = 
for all a E: A. 

If X, Y are random matrices, we wih let <p2 {X, y) = E (Tr (X) ,Tr{Y)). 



8 



Definition 2.14. Subalgebras Ai,. . . ,An of a second-order noncommuta- 
tive probability space {A,(fi,(p2) are complex second-order free if they are 
free and 

ip2 (ai • • • Up, 61 • • • 6p) = ^ V'l {a,ibk-i) 

fc=Oj=l 

when the ai , . . . , Op and the bi, . . . ,bp are centred and cychcally alternating, 
and 

ip2 (ai • • • Qp, 61 • • • 65) = 

when p ^ q and the ai , . . . , Op and the bi, . . . ,bq are centred and either 
cyclically alternating or consist of a single term. 

Here and elsewhere we take indices modulo the range over which they 
are defined. 

Definition 2.15. We say that random matrices {Xx : 0, — t- Mjsixn i^)}\(^A 
have a second-order limit distribution if there exists a second-order proba- 
bility space {A, (pi,(p2) with xx e A, X e A, such that for i = 1,2,..., any 
nonnegative integers ni,n2, ■ ■ ■ and any polynomials pi in noncommuting 
variables (i = 1, 2, . . .), and any Aj^i, . . . , Aj^„ € A, we have 

(^^ (^^I'l > ■ ■ ■ ' ^Ai,„i ) ) ) = (pi (a^Ai,! , • • • , a;Ai,„, ) ) , 



jJToo ('^ (^Am > ■ ■ • > J ) , IV (p2 [Xx,^, , . . . , Xa2,„, ) ) ) 

= V2 {pi (xAi,! , . . . , XAi,„ J , P2 (xA2,i , . . . , a;A2,„2 
and for all r > 3, 

iToo (^^ (^^I'l ' • • • ' ^"1-1 )),•■•> IV (Pr (Xa.,1 , ■ • . , ^A.,„. ) ) ) = 0. 

In addition, if there are A G A with —A G A, then we require the existence 
of an involution x 1-^ x* on ^4 reversing the order of multiplication such that 
x_A = x\ for any such A. 

Definition 2.16. For each colour c G /, let [x'^'^ : Q MnxN 

I J AeAc 

be an ensemble of random matrices. Let v : \p] ^ I and w : [q] —> I he words 

in the colours. Let A^ be a matrix in the algebra generated by the X^^^~^ 
and let -B^ be a matrix in the algebra generated by the X^^j^y Then we 
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say that the ensembles are asymptotically complex second-order free if they 
are free, have a second-order hmit distribution, and the algebras generated 
by the elements of the second-order limit distribution are second-order free; 
that is, for v and w cyclically alternating and p = q > 2, 



lirn^ k2 (tt (ii • • • ip) , Tr • • • Bp 



p-i p 



Vn 1™ {E{tT{ABk^i))-E{tT{Ai))E{tv{Bk_,))) 

^ — * A/ — 



fc=Oi=l 



and whenever p ^ q and v and w are either cyclically alternating or have 
length 1, 



lunj2 [Tr[A,...Ap),Tr[B,.--B,))=0. 

We collect the analogous definitions for real second-order probability in 
Section [71 



3 The matrix models 

In this section, we present formulas for the expected values of products of 
traces of matrices from the three models we are considering. Quantities 



involving several independent matrices can be calculated using Lemma 4.2 
We will use several lemmas to derive the formulas in this section, which 
we state here. The following is a folklore result which can be proven by 
direct calculation: 

Lemma 3.1. Let tt £ S {{ni, . . . , Um}), ni < ■ ■ ■ < Um < n, and let Xn^. be 
an Nfi^ X N^(^^^^ matrix for each k. Then 



Tr^ {Xn, , . . . , Xn^ ) - n ^iZl 

l<ini^<N„^^ k=l 
k=i,...,m 



If Nk = N for all /c G /, we can interpret this as a sum over all functions 
i : {ni, . . . ,nm} [N] with k i— )• ik- 

The following lemma is known as the Wick Formula. A proof can be 
found in |14], p. 164. 

Lemma 3.2 (Wick). Let {fx}^,^/^ be components of a real Gaussian family 
of random variables for some index set A. Then, for Ai, . . . , A„ € A, 

ie(/a,---/aJ= E n IE(/aJaJ. 

7reP2(n) {k,l}eTT 
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We will find the following expression useful. For e : [n] — >• {1, —1}, we 
define e on — [n] by e {—k) = e (k), and we define 



5s : k e (k) k. 

Lemma 3.3. The cycles of SsTrSnde are of the form (k, —s (k) e (l) I) and 
{—k,e(k)e{l)l), where {k,l) is a cycle ofn. It is thus a premap. 

Proof. The permutation 6 conjugated by vr consists of cycles of the form 
{k, —I), where {k,l) is a cycle of vr. Conjugating by 5s may interchange k 
and —k as well as / and — /. If £ {k) = £{l), then the cycle is unchanged; 
otherwise, the cycles become {k,l) and {—k,—l) (so no cycle contains both 
k and —k, and hence it is a premap). □ 



We will also use the following lemma, outlined in 20 . We give a complete 
proof here. 

Lemma 3.4. Let X be an M x N matrix with entries -4^ fa, where the 
fij are independent N (0, 1) random variables, and let Yi, . . .Y^ be random 
matrices independent from X and appropriately sized so the matrix multi- 
plication below is defined. (We will take the normalized trace to be times 
the usual trace regardless of whether it is applied to an N x N matrix or an 
M X M matrix; we will assume that M is of the order of N.) 
Let 7 G 5„ . Then 



= ^ ^#(7Z''5.vr57r5,7+)/2-#(7)-n/2 

Proof. To allow for the arbitrary occurrences of transposed copies of X 
encoded in e, in place of the usual indices ik and jk we will use indices z.^ 
and . We choose our indices so that the feth occurrence of the matrix X 
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will have indices t^i^ 



■ , ^^n V '■1 H 



t+:[n]^[Af] 
i-:[n]-5>[M] 

^ e{n) -e(n)-f -£(n) £(7(n)) 



(1) 

(1) ,£(7(1)) 
'7(1) 



^ \ £(7(1)) -e(n) e{7(")) I 

xE (/^+^- • • • f^+^- 



We now apply Lemma 3.2 to the expected value expression: 



vreP2(ri) {k,l}eTT 

If / or 7^ (that is, if t+ (/c) 7^ t+7r {k) or t~ (fc) 7^ l^tt (fc)) then 
E(/+-/+-) =0, so a product including this term will vanish. Thus the 
term corresponding to a given pairing vr is 

1, = o vr 



n ^ (■^^'-fc -^^^i ) {0, otherwise 



{k,i}eTT 

Inserting this back into our expression, we get: 



^ \ ^^-'(l)^e(7(l)) ^-e(n)^e:(7(n)) ^ ' 

,+ :[n]^[N] ^ ' " ^ 7re7'2(n):t±o,r=i± 



t-:[n]-s>[M] 

Reversing the order of summation, we get: 



i-:[n]-s>[Af] 



(1) ...V'W 

e(") e(7(n)) 
1 -7(1) S(n) 



By Lemma 3.3, SsTrdirSs and therefore 7_ (5e7r(57r5£7+ are premaps. We 
give an expression for the element of Yfc which appears, regardless of the 
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sign appearing on k in 7_ d^TTSTrde"/^: for k > 0, 

-sgn(fe)e(7_(fe)) sgn(fc)£(7_,_(fe)) -e(fc) e(7(fc)) 



|7-(fe)| |7+(fe)| 



7(fc) 



and 



_sgn(~fc)e(7_(-fc)) sgn(-fe)e(7+(-fe)) ^e(7('=) ) ^-e(fc) ^-e(fe) ^e(7(fe)) " 

Since exactly one of k and —k appears in the permutation 'yZ^SsTT6TTSsj+/2, 
we can write the expression: 

/ \ 



7rGP2 L+:[n]^[N] 
i-:[n]^[M] 



gn(fe)e(7_{fc)) sgn(fe)e(7^(fc)) 

\^A;e7l^(5E7r(57r(5E7+/2 "h-C^)! |7+(fc)| J 



We show that the condition o vr = always pairs a first index with 
a second index: specifically, the second index of Ifc is paired with the first 
index of Y^-i^ ^g^^ 7+(A:)' '^^lich must also appear, since 'jZ^SeTT67r6s'j+ (k) 

appears in the same (particular) cycle of ■jZ^Se'^^'^^el+- The second index 



3.3 



of Yu is L^^[^\'^S'^'^^''^\ which is paired with By Lemma 

|7+(fc)| ' J" ■^(n+{k)\) ^ 

JeVrJvr^e contains the cycle 

(7+ {k) , -sgn (7+ (k)) e (7+ {k)) e (I7+ {k)\)) ^ (I7+ {m , 

so TT (|7_|_ (A;)|) = \ j^jZ^ Ss7t6tt5s^+ {k)\ and (since the 7± do not change the 
sign) sgn {'jZ^6en6Tr5e^+ (k)) = -sgn{k) e {'y+ (k)) e {6eTr6TT5e'y+ (k)). Since 
each appears at most once in the expression, there are no further con- 
straints. 

Renaming the first indices and replacing the second indices with the 
index they are constrained to be equal to, Yk appears with indices ik and 
\-'^s.TTSiT5i:'y+{k)- recognize the resulting sum as a product over the cy- 



cles of 'j_^5s7rS7r6e^+/2 as in Lemma 3.1 and expressed in terms of the 



normalized trace, the result follows. □ 
Remark 3.5. It is possible to interpret the summands in these and subse- 



quent calculations as unoriented surfaced hypergraphs (see Remark 2.4). 
Briefly, the cycles of 7 can be thought of as encoding face information, and 
7+7Z^ face information on the two-sheeted covering space which allows for 
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a consistent orientation. The cycles of a premap vr encode hyperedge in- 
formation on this covering space. Vertex information on the covering space 
is given by the permutation 7+7r~^7Z^, and the traces of the Yk matrices 



are given by its inverse. See [14,16,20,25 for examples of how matrix inte 



grals may be interpreted in terms of surfaced hypergraphs. The definition 
of the Euler characteristic of 7 and vr which follows is the natural one in this 
interpretation ( \7\ gives a compatible definition of genus). 

Definition 3.6. Let / be a finite set of integers which does not contain both 
k and —k for any k. For a 7 G 5 (I) and a premap vr G PM (±/), we define 

X (7, vr) := # (7+7I') /2 + # (^) /2 + # (7;'^"S-) /2 - |/| • 

We note that if ±/i and ib/2 are disjoint, and 7^ G S (li) and vTj G 
PM{±Ii) for i = 1,2, then 

X(7i,7ri) + x{l2,TT2) = X(7i72,7ri7r2) . 
3.1 Real Ginibre matrices 

Definition 3.7. Let fij be independent A^(0,1) random variables, for 1 < 
hj ^ Let Z : i7 — )• M^xN (IK) be a matrix-valued random variable such 
that the ijih. entry of Z is fij . Then Z is a real Ginibre matrix (see ) . 

Lemma 3.8. If Z is a real Ginibre matrix, 7 G Sn, e : [n] — )• {1,-1}, and 
Yi, . . . ,Yn are random matrices independent from Z, then 



E tr 



Tre{p5p:pe'P2in)} 



Proof. By Lemma 3.3 dgTrde = depdpdg is a pairing, so # {dgTrSs) = n, and 



the lemma follows immediately from Lemma 3.4 □ 



3.2 Gaussian orthogonal matrices 

Definition 3.9. Let X : Q ^ M]\[xn (^) be a random matrix such that 
Xij = -^fij, where the fij are independent A^(0, 1) random variables. 
Then a Gaussian orthogonal matrix (or Gaussian orthogonal ensemble or 
GOE matrix) is a matrix T := ^ (X + X^). (This definition is equivalent, 
up to normalization, to more standard definitions given in, e.g., [l5|. See 
e.g. [8, 14 for a demonstration of the equivalence in the complex case; the 



real case is similar. 
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Lemma 3.10. If T is a GOE matrix, 7 G Sn, and Yi, ... ,1^ are random 
matrices independent from T , then 

E(tr^ {TYi,...,TYn)) 

7rePM(±[n])n7'2(±H) 

Proof. We expand the left-hand side in terms of X and and apply 
Lemma 13.41 to each summand: 

E(ti>(ryi,-- - ,ry„)) 

= _L V V /V#(7l''5.P'5p'5.7+)/2-#(7+7l')/2-n/2 

2n/2 

£:[n]^{l -l}peT'2(n) 

xE ftr -1 



We now show that the map (/), e) 1— >■ d^pdpSe is 2"/^-to-one and its image 
is PM (lb [n])n'P2 (=t [f^])- By Lemma [3.3| SsP^pSg is a pairing and a premap. 
For any vr € PM (it [n]) n P2 (i [n]), we can find a p £ V2 (n) and an 
appropriate e such that we construct vr in this manner. In fact, we have two 
choices for values of e on the elements of each cycle of vr, or 2"/^ such e. 

Since vr is a pairing, # (vr) = n, giving us the correct exponent on N, 
and the result follows. □ 

3.3 Real Wishart matrices 

Definition 3.11. Let {fij} be independent N (0, 1) random variables, for 
1 < i < M and 1 < j < iV. Let X : MmxN {^) be a matrix- 
valued random variable such that the ijth entry of X is ~^fij- Then a real 

Wishart matrix is a matrix W = X^DX for some M X M matrix D. We 
will assume that M is of the order of N, and each D is part of a family 
Df^^ G MAfxAf (C), A G A such that limiv^oo tr (^D^f ^ • • • exists for 

all n and all tuples Ai, . . . , A.„. 



Remark 3.12. Wishart matrices were first studied in 24 , and are variously 
defined in such papers as [l|[3|[4|[lT|[l2|[T6]. They are generally defined as 
matrices of the form AX^ BXC with various requirements on the constant 
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matrices A, B, and C, such as that they must be symmetric, positive def- 
inite, or the identity matrix, or various relationships between them. Here 
we require A and C to be the identity matrix, but we do not put any re- 
quirements on B. The possible asymmetry of B and therefore the Wishart 
matrices motivates the appearances of the linear involution in the main def- 



inition 7.2, and the possibility of different B matrices motivates the order 
in which terms appear. 

The following lemma expresses the moments of real Wishart matrices in 
a form which we will find convenient. The pairings on [2n] which appear in 



Lemma 3.4 are replaced by premaps on it [n] under a bijection in which each 
pair of adjacent elements is compressed into one element. This is possible 
because the transposes appear in a regular pattern; in some sense, sign takes 
the place of parity in this expression. Geometrically, this calculation can be 
represented as in [20] , except we consider the alternate vertices which here 
contain the matrices as hyperedges, encoded in the permutation vr, and 
the vertices induced in this unoriented hypermap are the vertices containing 
the Yfc matrices. An example with a diagram is given after the proof. 

Lemma 3.13. If{Wx : 0, — t- M]\fxN {'^)}xeA ^'^^ '^^^^ Wishart matrices such 
that Wx^ = X^Dx^X for 1 < k < n, j £ Sn and Fi, . . . , y„ are random 
matrices independent from the Wxj, , then 

K{ir.,{Wx,Y^r-- ,WxYn)) 

= ^ m(^'-^-'n^hT^_,^Dx,,...,Dxj 

7T£PM{±[n]) 



7- 'r7+/2 

Proof. In order to put the left-hand-side of our expression into the form of 



Lemma 3.4, we define two bijections from the indices [n] to the new indices 
of the X-^ and X terms (respectively) that appear when each Wishart term 
is expanded: 

6- :k^2k-l. 

and 

9^:k^2k 



We then define, for k G [2n] 
I'ik)-- 



e+ez^ (k) k odd 

e--f97\k) A; even ' 
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e:k^ {-ir, and 5' = de- Then 

E(tr^(VFA,yi,-- - ,WxYn)) 
= E (try [X^D^.^XYi,..., X^Dx^ , Xy„) ) 

The permutation j'S^S' pSp6'j'_^_ acts separately on the odd and the even 



integers: by Lemma 3.3, 6'p6p5' reverses either the parity or the sign. In 
the former case, exactly one of 7^ and j'S^ acts, and in the latter both or 
neither do. In either case, the parity is preserved. 

We extend the 6± as odd functions (so they commute with 5, 6', and 
absolute values). We define vr G PM (± [n]) by vr"^ = OZ^'y'S^S' p6p6'-/'^9-. 
Then the contribution of the odd cycles is tr^-i /2 (-^An ■ ■ ■ , D\^). We now 

show that the contribution of the even cycles is E {^'^^-'^^^^^12 (^i) ■ ■ ■ t Yn)^ ■ 

We note that 7'^ = 0±'^0^ for whichever of 9^ is defined, and O+OZ^ = 
7'(_7'_. Then 

^+71^7+^;^ = j'S^9+9z\+^6' p6p6'-f'_e-e^^j'^ = ilH' p5p5'i+. 

We can see that 7^ (7') = 7^ (7), since odd integers always occur in cycles 
with even integers, and on even integers, 7'^ = ^+7^^^, so the exponent on 
N is correct. 

We now show that this map from V2 (2n) to PM (it [n]) is a bijection. 
For a given vr G PM (it [n]), we can reconstruct from vr and 72^7+ the 
permutation j'S^S' pdp6'^'_^^ and hence 6'p6p5', and by Lemma 3.3 we can 
reconstruct p. Thus the map is injective. 

Conversely, for vr G PM (it [n]), we can define /3 on /c G [2n] by letting 

P{0± (k)) = |6'=Fsgn(7rTi{fc))^'^^ ' 

(All integers in [n] are in the range of exactly one of the ^±1.) If p {0± (k)) = 
d±i (k), then sgn (vr^^ (/c)) = —1 and Ivr^"*^ {k)\ = k, which is excluded since 
TT is a premap; thus p has no fixed points. We calculate: 

p' (0±i {k)) 

^±sgn(^Tl(fc))sgn(7r±=S"(-^'W)(|7rTl(fc)|))^^'^''^^^ ^^^^ d^^^ D 
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Figure 2: A pairing (left) and its equivalent premap (right). "Twists" in 
the cycles of the premap, which denote that that edge is identified in the 
opposite direction, correspond to changes in sign. 



Compensating for the absolute value signs, we note that 

^±sgn{.^Hk)) (I^Tl = ^±sgn(.Ti(fc)) ^^^^ (^^1 (fc)) ^Tl (fc)) 



Sg 



sgn ITT 



(k)) k, 



we can see that p is the identity. Thus p is a pairing. 

Finally, we show that 0Z 7^ 5' p5 p5' ^'__9 - is premap vr. Since the permu- 
tation 'y'^^ 5' p5 p6' takes odd numbers to odd numbers, ^'_9- = 6'_sgn(jfc) 

-' = 0, 

sgn{k) (k) I 



and 01 



0^^ for whichever 9^^ is defined. We have 



sgn(A;) 



sgn ( A: ) sgn (tt ( ) { I A: I ) ) 

^sgn{A:)2sgn(7r(fe))^ (^) 



^sg 



n(7r{fc))7r (^)| 



SO by Lemma 3.3 6' p5p5' takes 0— sgn(A;) 

(A;) to 

- sgn (6'_sgn(fc) (A;)) e {9_sgn{k) (^)) ^ {P (|6'-sgn(fc) (^)|)) P (l^-sgn(fc) (^)| 

= -sgn (k) e {9-sgn{k) {k)) e (6'sgn(7r(fe))7r {k)) \9sgn(n(k))'^ (k) \ 



= sgn(/c)^sgn(7r(A;)) |6'sgn(7r(fc))^ (^)| 
The result follows. 

Example 3.14. If we wish to calculate the value of 



'sg: 



:n(fc)sgn(7r{fc))7r 



□ 



E (tr iWiYiW2Y2W3Y3) tr (WiYiW^Y^)) 
= E (tr {X'^ DiXYiX'^ D2XY2X'^ D3XY3) tr {X'^ D^XY^X'^ D5XY5)) , 
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Lemma |3.4| gives us the expression 

P6P2([10]) 

where 

7 = (1,2,3,4,5,6) (7,8,9,10). 

Each term in the sum can be represented as a face gluing (see [20], Section 3, 
for more detail). The pairing 

/9= (1,5) (2,7) (3,9) (4,10) (6,8) 

is shown in Figure [2| left. We calculate that 

(1, 7, -5) (5, -7, -1) (2, -8, 6, -4, 10) (-10, 4, -6, 8, -2) (3, -9) (9, -3) . 

The region that becomes the vertex (1, 7, — 5) (5, — 7, — 1), which becomes 
tr (^DiD^DJ^ , is marked with arrows. The contribution of this term is then 

N-hi [DiD^Dl] tr {D2DI) E (tr (YiYIy^Y^Y^)) . 

As shown in Figure 

§ right, we can "pinch" the two edges corresponding 
to each Wk matrix together. We now have 

7 = (1,2,3) (4,5), 

and Lemma |3 . 1 3| gives us that 

E (tr {WiYiW2Y2W^Y:,) tr {WiY^W^Y^)) 

= ^ mM-hT^-./^ iD,,...,D,)E (tr^-1,,^/2 (Yi,..., n)) 

7rePM(±[5]) 

The above pairing p corresponds to premap 

^ = (1,-3, 4) (-4, 3,-1) (2, -5) (5,-2). 

The cycles (1, — 3, 4) (— 4, 3, — 1) are the hyperedge corresponding to the 
marked vertex in the previous diagram. This premap has vertex permu- 
tation 

7li^7+ = (1, -4, 3, -2, 5) (-5, 2, -3, 4, -1) . 
We calculate that 

X (7, vr) = 2 + 2 + 1-5 = 0, 
so vr gives the same contribution as p did above. 
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4 Combinatorial calculations 



Let {-^^aIasa ^^^^ Ginibre, GOE, or real Wishart matrices. (We consider 
a family of matrices in order to accommodate Wishart matrices Waj, = 
X'^Dx^X with the same matrix X but possibly distinct deterministic ma- 
trices D\^\ in the Ginibre and GOE cases there is only one matrix Xx. 



Independent ensembles will be handled by Lemma 4.2 ) We note that for 
all n, 7 G S*,!, e : [ra] — {1, — 1}, and random matrices independent from 
the X\, the Xx satisfy 



E^tr, [Xl'^'^Yu--- ,Xll''»Y„ 
= ^ iV^(^.^-^^)-2#(7)/^(,,)E(tr^-.,^^,^^^/2(yi,...,y„)) (2) 

7rePM,(±[n]) 

where PMc{i[n]) is a subset of the premaps on ± [n] and fc is a func- 
tion on those premaps. In each case, for each finite set of positive in- 
tegers /, PMc (H) ^ PM (±/) is a subset of the premaps on ±7 such 
that for any J Q I, the vr G PMc (H) which do not connect it J and 
± (/ \ J) are the product of a vri G PMc (±J) and it2 G PMc (± (/ \ J)), 
and fc ■ U/cN |/|<oo (i^) — 5- C is a function such that lim7v-!>oo fc (tt) 
exists. Furthermore, if vr G PMc{I) does not connect it J and ib(/\ J), 

then fc (vr) = /c {'t^\±j) fc (^^l±(/\j))- (This last condition is the only one 
not satisfied by Haar-distributed orthogonal matrices (appearing in future 
work). As it is only required for some of the results, we will note when it is 
needed.) 

Specifically, for real Ginibre matrices, PMc (±-^) = {p^p : p G 7^2 {!)} 
and fc (vr) = 1. For GOE matrices, PMc (±7) = PM (±7) n V2 (±7) and 
fc (vr) = 1. For real Wishart matrices, PMc (±7) = PM (±7) and, if Xx^ = 
Wxk = X'^Dx^X and 7 = {n, . . . with ii < . . . < im, then fc (vr) = 

tr^-1/2 {D\,^ 7)a,„ 



4.1 Moments and cumulants 

We now give a formula for the traces of products of matrices from several 
independent ensembles of matrices satisfying ([2]). It does not depend on 
the last condition (that / be multiplicative). This formula can be used to 
calculate expressions with several independent Ginibre, GOE, or Wishart 
matrices, or any combination of such matrices. 
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Remark 4.1. Lemma |4.2| is the fairly intuitive result that, given a product 
of traces of several independent ensembles of random matrices, we construct 
faces as before and glue the edges belonging to each independent ensemble 
according to the rules of that particular ensemble. An example with a 
diagram is given after the proof. 

Lemma 4.2. For each colour c G [C], let \Xc^^\ be a set of random 

matrices satisfying ^ with subsets of the premaps PMc {^I) and function 
fc ■ U/cN |/|<oo -^-^c (i-f) C, and assume the set associated with each 
colour c is independent from every other set. Let w : [n] ^ [C] be a word in 
the set of colours [C]. Then, for 'j £ Sn and e : [n] — t- {1, —1}, 

E(tr,(x«;^^-),....xitr»')) 

iV'<'^*"'->-'*<''/i(>ri)---/c(»o). (3) 

7rcePMc(±«)-l(c)) 

Proof. We prove this lemma by induction on the number of colours C. For 
C = 1, Formula p]) reduces to where Yi = ■ ■ ■ = Yn = In- We now 
assume Formula pf for C — 1 colours. 

The left hand side of ([s]) is of the form of Q in colour C, where 
7' := 7|^-i((7) takes the place of 7. (We may cycle any matrices appearing 
in a trace before the first occurrence of a matrix of colour C to the end with- 
out changing the value, and any traces which do not contain any matrices of 
colour C may be treated as a constant multiplicative factor within the expec- 
tation.) Let w ^ (C) = {fci, . . . , kiuj- MC)|}' h< ■■■ < fc|«,-i(c)|- We can y^,^ 

the product of the matrices appearing (cyclically) between X^^ '''■^ and 



(e(fcy(i))Afc 

X^, ^ . (This may be no matrices, in which case we let Yk- := In-) 

For any trace which does not contain any matrices of colour C, we call the 
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matrix in the trace Zfc, 1 < A; < # (7) — # (7'). Then: 
xtr(Zi)---tr (Z#(^)_#(y))) 

^ JV^{7',<5e7rc5e)-2#(y) (^^) 

7rcePMo(±«;-i(C)) 

(tV-^4-c5.7;/2 (^A^i' ■ ■ • ' ^k^^-Hc)\) (^1) • • • (^#(7)-#(7'))) • 

We note that each matrix of one of the colours in [C — 1] appears exactly 
once inside the expected value expression, possibly transposed. Let us de- 
note by / the set of signed integers representing the indices of the ma- 
trices which appear. By the induction hypothesis, for some permutation 
7" G S {I), counting cycles as permutations on ±w~^ {[C — 1]), 

^ K-l5..c47;/2 (^fel' ■ ■ -'^^-1(^)1) ^""^^'^ (%(7)-#(7'))) 

^ JVX(V',4.'4)-2#(T")/^ (^^) . . . f^_^ (^^_^) . 

7r'=7ri-"7r(7_i 
7rcePMc(±«)-i(c)) 

(We note that PMc [±w~^ (c)) is still the appropriate subset of the premaps: 

if —A; G / for k > 0, then —k would appear instead of k in the premaps in 
PMc [±w~^ (c)), but the sign of e (k) and hence 6^ (k) would be reversed, 
so the permutation ^eVTc^e is unchanged.) 

Substituting this expression into the previous one, it only remains to 
confirm that the exponent on A'^ is correct. To do so, we examine the per- 
mutation 7". Wc show that 7" = 6siTcSel+lZ^\j- 

For i in a cycle of 7 which does not contain colour C, 7" (i) = 7 (i) = 

If a term with index i appears in signed, and hence i of the same 

sign as k), then the next term, if there is one, has index j+jZ'^ (i) (the 
correct permutation acts given the sign of i and hence the order in which 
the terms of appear). If there is such a next term, w (7+7!^ («)) 7^ C, so 
S^TTcSe acts trivially on 7+7! (i), and hence 7" (i) = SsTTc^el+lZ^l]- (i)- 

On the other hand, if there is no such next term, then 7" (i) is the first 
term of ^(^f'-^SeT^cSe'y' )'"(fc)' ^^^^^ smallest positive integer such 
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that this term is nontrivial. If i is the last index appearing in Yj-, then 
k = 7^17^7+ 7!^ (i) il'^^ acts depending on the sign of k, i.e. on whether k 
is to be found before or after the indices appearing in Y^). The first term 
of ^7l-i5,7rc<5e7;(fc) exists) has index j+-fZ^i_-f'_r^6eTTcSe7+7Z'^ (i) = 

7+7Z^(^£7rc(^£7+7Z^ (i)- If, however, this index has colour C, it is nonetheless 
the index encountered after the last term of Y,/-ic r ,/ in which case 

we may apply permutation j^jZ^ S^iTcSe until we have an index not of colour 
C. Since 7^ C exactly when w {\Sg-JTcSs (/)|) / C (and then deTrc^s 

acts trivially), the index of the next term is indeed SsTTc^el+lZ^lj {">■)■ 

On -/, -/'!r^ = 65eTTc6£j+jZ^5\_] = '^+'yZ^6£7rcSe\_j- Since SsTVcSs 
acts trivially on ±w~^ {[C — 1]), this is equal to 7"^"*^ = 6s7Tc^^+^Z^\_j- 
Thus 7+7""^ = deTTc^e^+lZ^l^p SO 5£7r'47"7""^ = 6sn5e^+^Z^\^j. 

We now show that the cycles of SeTrSsj+jZ^ which do not appear in the 
permutation it induces on ±7 = :izw~^ ([C — 1]), that is, those which are en- 
tirely colour C, arc in one-to-one correspondence with cycles of 7^ 6£Tt6£^'_^_ 
consisting entirely of k with trivial (that is, those for which it;(|A:|) = 
(7 (I^D) = C), which we count. We note that 7 (|A:|) = |7+7- {k)\. 

If a cycle of SeirSsj-^-jZ^ consists entirely of k with (|A;|) = C, then it 
is the factor ttc of vr which acts, so w (|7+7r^ (^)|) ~ ^- Since only one 
of 7+ and jZ^ acts nontrivially, w (I74- (A;)|) = w {{jZ^ (^)|) — so 
7+ (k) = 7+ (k) and j'S'^ (k) = 7_ (k). Thus this cycle also appears as a 

cycle of SeTTcSel'+l-^ ■ 

If we map each integer in the cycles of 6sTrcSsj'_^_j'S^ under 7^7^, they 
become the cycles of 7^^(5£7rc(5e7+. If /c is in a cycle of SeTTSej+j- which 
is entirely of colour C, then w ik)\) = C and w (|7+7_7^^ (/c)|) = 

w (I7+ {k)\) = C. Thus ^'ZZ^ {k) in the corresponding cycle of Seirc^el'+l-^ 
has Xy'-i(jfc) trivial. 

Conversely, cycles of ^'z'^^e'^C^el'+ correspond to cycles oi Sf^TTc^el'+l'-^ 
with each integer mapped under If Y^ is trivial, that is, if tt;(|fe|) = 
w (I7+7- {k)\) = C, then w {\-/'_ {k)\) = w (|7_ {k)\) = C. Thus a cycle in 
7^r^(5£7rc5£7+ of k with trivial corresponds to a cycle of (5£7rc(^£7+7^^ of 
colour C, which is equal to a cycle of SsTtSsJ+jZ^ ■ 

There are # (7) — 7^ (7') cycles of 7" which correspond to traces of 
Zfc matrices and hence # (7") -|- # (7') — # (7) which contain the Yfc and 
hence correspond to cycles of j'S^SeTTcSeJ+f^- Thus # (7Z^(5£7r(5£7+) /2 = 

# (,5£7r(5£7+7l'L^-i([c-i])) /2 + # (7- '<^.7rc5£7;) /2 - # (7") - # (Y) + 

# (7). The result follows. □ 
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Figure 3: The faces associated the calculation in Example 4.3 involving sev- 
eral independent ensembles and an example of an edge-identification con- 
tributing to the result. 



Example 4.3. Let C = 3. Associate with colour 1 the set of Wishart matrices 



where Xi is a random matrix as in Definition 3.11; with colour 2 an inde- 
pendent set of Wishart matrices 

[W^^^ = X^BiX2, = X^B2X2, ...} 

where X2 is another random matrix as in Definition 3.11 independent from 
Xi] and with colour 3 the Ginibre matrix Z, independent from the Wishart 
matrices. 

If we wish to calculate the quantity 

E (tr ( ZW^^'^) tr ( V^f ^^Z^Z^) tr ( H^f ^^Z^W^f « Vf 



for some integers A2, A3, Xq, As and Ag, then we let 

7 = (1,2) (3, 4, 5) (6, 7, 8, 9), 

we let w{3) = w (9) = I, w{2) = w (6) = u> (8) = 2, and w (1) = w (4) = 
w{5)=w (7) = 3, and we let e (1) = e (2) = e (3) = e (6) = e (8) = e (9) = 1 
and e (4) = e (5) = e (7) = —1. We construct faces shown in Figure [sj 
We then consider the hyperedges we can construct on edges 3 and 9, the 
hyperedges we can construct on edges 2, 6 and 8, and the ways in which 
edges 1, 4, 5 and 7 may be identified pairwise in the directions shown. We 
show an example of a choice for each corresponding to 



^1 = (3) (-3) (9) (-9) 
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(contributing a factor of tr (Ax^) tr (^Ag)) 



^2 = (2, 8, -6) (6,-8,-2) 



(contributing a factor of tr (^B\,^Bj^^B\g^), and 



7r3 = (1,-7) (-1,7) (4, -5) (-4, 5) 



or 



5,7r3(5e = (1,7) (-1,-7) (4, -5) (-4, 5). 



This gives us 



= (1, 7) (-1, -7) (2, 8, -6) (6, -8, -2) (3) (-3) (4, -5) (-4, 5) (9) (-9) . 



lZ^5eT^5el+ = (1, 8, 9, -7, -2, 6) (-6, 2, 7, -9, -8, -1) (3, -4, 5) (-5, 4, -3) . 



We now show that, if the last condition of Q (the multiplicativity of /) 
is satisfied, the cumulants of traces of products of matrices are sums over the 
premaps which connect blocks ±1, where / is an orbit of the permutation 
7. Such terms can be thought of as the connected surfaces. 



random matrices satisfying ^ with subsets of the premaps PMc and func- 
tion fc, and assume that the set associated with colour c is independent 
from each other set. Let w : [n] ^ [C] be a word in the set of colours 
[C]. Let ni,...,nr be positive integers, let n := ni + ■■■ + Ur, and let 
-^fc = ['^i + ■ ■ ■ + n-k-i + 1, ni + • • • + Tifc], 1 < k < r. Let 

7 = (1, . . . ,ni) • • • (ni H h Ur-i, . . . ,n) 

and e : [n] — t- {1,-1}. For 1 < k < r, define classical random variable 



We calculate that 



The contribution of this term is then 





be an ensemble of 



Yk ■= tr X 



.(e(niH hnfc_i +l)A„j 

ui(niH f-nfe_i+l) 



■■■X 



.(e(niH l-nfe)A„j+...+„j.) 

w(niH hrifc) 



) 
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Then 



7r^ePMc(±«i-l(c)) 
7rV{±4}^^i=l±[„] 

where the last condition under the summation sign means that vr must con- 
nect the blocks ±/i, . . . , zt/^. 

Proof. We demonstrate this by showing that the moment-cumulant formula 
([T]) is satisfied by the conjectured expressions. 

Consider the term corresponding to a vr in the moment 

7r=7ri ■■■7r(;7 

Any such vr induces a partition p €z V (r): we let p be the smallest partition 
such that s,t £ [r] are in the same block if zizlg and it/t are in the same block 
of TT V {ib/fc}^^-|^. Letting ly = Ufcey for V CI [r], the term corresponding 
to vr can be expressed: 

iV«<'*"'-'-'7l('ri)---/c(>rc) 

.niV«("'-'-"-''.')-l''l/.(..|,^.). 

yep 

We note that the induced permutations on ly are in each case simply 
the restrictions, since ±/y is the union of orbits in each case. In fact, 
for vr inducing partition p, the 7r|j^ are exactly the vr' = 7r^---7r^, tt^ G 
PMc {±w-^ (c) n ±/y) which connect the ±4, A; G F. The sum over all tt 
inducing p is then 

yep ,r'=7ri---7r^ 

<ePMc(±'LO-i(c)nJ'v) 
7r'V{±4:fcey}=l/^ 

where the last condition under the summation sign means that vr' must 
connect the ±1^ for all kinV. If ^ = {ii, . . . ,is}, then each term in 
the product is the conjectured expression for the cumulant kg (Yij, . . . , FjJ. 
This demonstrates that the conjectured expression for cumulants satisfies 
the moment-cumulant formula. □ 



26 



We will also find it useful to consider expressions in which terms in the 
algebras generated by each ensemble have been centred. A cumulant of 
products of centred terms may be interpreted in terms of the Principle of 
Inclusion and Exclusion, a generalization to functions of the formulas for 
calculating the number of elements in certain subsets and not in others (see, 
e.g. js], Chapter 5, for more details and a proof). 

Lemma 4.5 (The Principle of Inclusion and Exclusion). Let 



This allows us to write a cumulant of products of centred terms sat- 
isfying the last condition of ^ (multiplicativity of /) as a sum over the 
permutations with certain connectedness properties. In addition to being a 
connected surface, each interval of edges corresponding to a centred term 
must be connected to at least one other. 

Lemma 4.6. Let pi, . . . ,Pr be positive integers and let p := pi + ■ ■ ■ + pr- 

Let w :[p\ ^ [C] he a word in the set of colours [C]. Associate with each 



Let ni, . . . , np he positive integers, let n := ni + ■ ■ ■ + Up, and let 
7 = (1, . . . ,ni H h npi) (ni H h n-p^ + 1, . . . , ni H hnpj • • • 



f,g:{IC[n]}^C. 



Then the following are equivalent: 



g{i) = Y,f{J); 



JDI 



/(/) = E(-i)"^"5(^)- 



colour an independent matrix ensemble 





let e : [n] — )■ {1, —1}, and let 



niH hrifc 



n 



X 



-(e(i)Xi) 
w{k) 



i=ni-\ hrtfc-i+1 



For k £ \p], let 



Ik ■=[ni-\ h Uk-i + 1, ni ^ h nfc] . 
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For K C [p], let Ik = Ufcex^fc- Vk = /[p,+...+p^._j+i,pi+...+p^] . Then 



fcr ( tr ( Ai • • • J , . . . , tr ( Ap,+...+p^_i+i ■■■Ap^ 

J2 ^x(7A-4)-2r^^(^^)...J^(^^) 



7rV{±ya[^l = ln 
±/fe^7rV{±/fc}^^l 

(^f/ie t/iird line under the summation sign says that n connects the blocks 
ibVfc, and the fourth that it connects each block zt/^ to at least one other). 

Proof. Expanding the left-hand side expression, we get 
(-l)l^lniE(tr(A,)) 

KC\p] keK 

xkrltrl Y[ AJ,...,tr Y[ Ak 

\ \k(i[l,pi]\K 1 \fee[pi+-+Pr-l + l,p]\X 

For each A; G [p], define the single cycle permutation 

7A: := (ni H h nfc_i + 1, . . . , H h n^) . 

Applying Lemmas |4.2| and |4.4| to our previous expression, we get 



^ — ^ )('( 7lr I'^e'T'Je |— 2r 



7rV{±yi}[^l = l±[„j 

(We note that if any [pi + • • • + Pfc-i + 1, pi + • • • + Pfc] \ is empty then 
one of the entries of the cumulant is tr (/at) = 1, so the cumulant is zero and 
the term can be ignored. Thus we can assume that # (^l/f y\k) ~ ^'-^ 

For a given the permutations in each sum act on disjoint sets, so we 
can express it as a sum over permutations on all of ± [n] . Let 



7' 
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Then 



_ft:C[p] keK 7r=7ri---7rc 

7rV{±ya[^l = l±[„j 
±7fce7rV{±7fc:fce/^} 

(where the last condition under the summation sign means that tt does not 
connect ±7^ to any other ±7/, for all k £ K). 

We now show that the exponent on is equal to the exponent if permu- 
tation 7' were replaced with 7. We recall that multiplying a permutation tt 
by a transposition {k, I) (on the left or right) joins the orbits of tt containing 
k and I if they are in separate orbits, reducing the number of orbits by 1, and 
splits the orbit containing k and / if they are in the same orbit, increasing 
the number of orbits by 1. 

We note that 

7 = 7 n ("i + --- + '^fc'7"^|[„]\;^("'i + --- + rifc)) 

(straightforward calculation; the product of transpositions is in increasing 
order in K and inverting the product reverses the order of the transposi- 
tions). Each transposition connects an originally disconnected interval to 
an orbit of 7', so # (7) = # (7') — \K\. 

We then express ^Z^5eTT5e^+ as ^'z'^SeTide'^'j^ left- and right-multiplied 
by transpositions. The permutation 7^^5e7r(5e7^ does not connect any 1^ or 
—Ik for G -fC to any other set, so each of these transpositions reduces the 
number of cycles by 1. Thus # (7r"^(5e7r(5e7+) = # {^'Z^^eT^^e^'-^.) — 2 \K\. 

Changing the exponent on N accordingly, we have: 

KQ\p] 7r=7rr--7rc 

7rV{±ya[^l = l±[„j 
±/fee7rV{±/fe:fce/r} 



We can then interpret this expression in terms of Lemma 4.5 In our 
case, we will let /, (7 : {K C [p\] — )■ C, with g [K) the sum of terms over tt 
which do not connect the intervals ±1^ to any other interval for all A; G 
and / {K) the sum of terms over vr which do not connect the intervals ±1^ to 
any other interval for exactly the k £ K (i.e. any other interval is connected 



to another). We note that / and g satisfy the hypotheses of Lemma 4.5 



The desired quantity is equal to / (0), from which the result follows. □ 
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5 Asymptotic calculations 



We may find upper bounds on the order of N and characterize highest-order 
terms in many of the above formulas using the following well-known result 
(see, e.g. (lO|, and [l6|] for a proof): 

Theorem 5.1. Let n, p £ S (/) for some finite set I. Then 
# (vr) + # (vrp) + # (p) < |/|+2#(^,p). 
We see that if 7 G and vr G PM (it [n]) represent a connected surface 



(see Remarks 2.4 and 3.5), then x (Tj ^) ^ 2, as we would expect of an Euler 
characteristic: 

Lemma 5.2. Let 7 G Sn, and let {Vi, . . . , Vr} G V (n) be the orbits of If 
IT G PM (lb [n]) connects the blocks of {±Vi, . . . , ±Vr}, then x {l-, t^) ^ 2. 

Proof. If an orbit of (7+73^, tt) does not contain at least one of or —V^, 
then vr does not connect the block ibVfc. Thus 7^(7+73"^, vr) < 2. The result 
follows. □ 

5.1 Highest order terms and noncrossing conditions 

For a given vr, a p satisfying the equality (typically the ones which will 
contribute highest order terms in N) may often be interpreted as a non- 
crossing diagram on the cycles of vr. For vr with one or two cycles and (vr, p) 
transitive, we state conditions equivalent to satisfying the equality. See the 
original references for proofs and diagrams. 

Remark 5.3. We are following the conventions of [7, 14 , where vr is thought 



of as enumerating the hyperedges, which should also be counterclockwise; as 
opposed to those of p]p[6], where a standard vr should follow the sense of 7. 
The vr which here satisfy or violate the various conditions are the inverses 
of the vr of the latter sources. 

Definition 5.4. Let 7 G S {I) for some finite set I be a permutation with 
a single cycle, and let vr G 5" (/). 

We call vr disc nonstandard (relative to ^) if there are three distinct 
elements a,b,c £ I such that 'jl^abc} ~ (^!^)C) and ttI^^ be} ~ (^)^)C). We 
call vr disc standard (relative to 'j) if there are no such elements. 

We call vr disc crossing (relative to 7^ if there are four distinct elements 
a,b,c,d e I such that 'y\^a,b,c4} = ("■■• ^' ^) ^ut Tr\^a,b,c,d} = '^)- 
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call vr disc noncrossing (relative to if it is neither disc nonstandard nor 
disc crossing. 

We denote the set of disc-noncrossing permutations on / relative to 7 by 

'S'disc— nc (t) ■ 

The following theorem is from (2|: 

Theorem 5.5 (Biane). Let ^,it £ S (I) for some finite set I, where 7 has 
a single cycle. Then 7 and vr satisfy the equality in Theorem 5.1, or in this 
case 

#(7r) + #(7-V-i) = |/| + 1, 

if and only if n £ S'disc-nc (t)- 

Similar conditions may be found for permutations connecting two cycles: 

Definition 5.6. Let 7 G S (!) for some finite set / be a permutation with 
two cycles (which we will refer to as 7ext and 7int), and let vr G 5 (I). 

We say that vr is annular nonstandard (relative to 'j) if one of the two 
following conditions holds: 

1. there are a,b,c £ I such that 7||„ be} ~ i^^ ^' ^^'^ ^l{a be] ~ ('^' ^1 '^)' 

2. there are a,b, c, d £ I such that "y\^abcd} ~ C*^' ^) i'^^ d) but 7r| |^ ^ ^ ^| = 
(a, c, b, d). 

We call TT annular standard ( relative to 'j) if neither of these conditions holds. 

Let X £ 7ext and y £ 7int- We define a permutation X^^y on I\ {x, y} by 
letting X^^y (7"^ (x)) = 7(2/) and Xx,y {j^^ (y)) = 7 (x) (we will generally 
be assuming that 7ext and 7int each have at least two elements) , and letting 
'^x,y (o) = 7(0) otherwise. We will say that vr is annular crossing (relative 
to 'J ) if one of the three following conditions holds: 

1. there are elements a,b,c,d £ I such that jl^abcd} ~ {<^-,b,c^d) but 
7r|{a,b,c,d} = (6,d), 

2. there are elements a,b,c,x,y £ I, x £ 7ext and y £ 7int, such that 
^^,y\{a,b,c} = («>^c) and = ia,b,c) {x,y), 

3. there are elements a,b,c,d,x,y £ I, x £ 7cxt and y £ 7int, such that 
>'^,y\{a,b,c,d} = {a,b,c,d) but vr||„ = (a,c)(6,d) (x,y). 

We call vr annular noncrossing ( relative to 'j) if it is neither annular non- 
standard nor annular crossing. 

We denote the set of connected annular-noncrossing permutations (those 
for which (7,7r) is transitive) by 5ann-nc (7)- 
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The following theorem is from p!6|: 



Theorem 5.7 (Mingo and Nica). Let 7 G S {I) for some finite set I he a 
permutation with two cycles, and let tt & S (/) with it connecting the cycles 



of J (i.e. #(7,7r) = 1). Then 7 and it satisfy the equality in Theorem 5.1 
or in this case 

#(vr) + #(7"V-i) = |I|, 
if and only if n £ S'ann-nc (7) • 

We generalize these results to the nonorientable case. 

Lemma 5.8. Let 7 G 5 (/), where 7 has a single cycle. Let vr G PM (±/). 

Then X (7) ^) = 2 i/ and only if tt does not connect I and —I and ■n\i G 

'S'disc— nc (7) • 



Proof. If x(7)'^) = 2, then by Theorem 5.1, #(7+7_^,7r) > 2 so vr must 
not connect the orbits of 7+7!^, / and — /. Thus 7+"'^vr^^7_ j^^ = ^^^-k^'^, 
and since vr and 7J|I"^7r~^7_ have the same number of cycles on / as on —I, 
+ # (7"'vr-i|^) = |/| + 1. Thus G 5disc-nc (7)- 
Conversely, if vr does not connect / and — / and vrlj- G S'disc-nc (7)) we 
calculate similarly that x (7) ^r) = 2. □ 

Lemma 5.9. Let 7 G S{I), where 7 has two orbits, V\ and Vi- Let vr G 
PM (±/) connect itVi and ibV2. Then x (7, tt) = 2 if and only if, for some 
choice of sign e = ±1, vr does not connect Vi U eV2 to (— Vi) U {—€¥2) and 

■^IviUeVa ^ •S'ann-nc ( 7+7- | y^ueVa, 

Proof. Assume x (7i^) = 2 and tt connects itVi and ibV2. Being a premap, 
the permutation vr must connect Vi to either V2 or — V2 and —Vi to the 



other. We can calculate from xili'^) s-^id Theorem 5.1, however, that 
7^(7+7-^) — 2, so vr, and hence 7^^vr^^7_, must not further connect these 

blocks. We calculate that # \ AvMeV2)) ( ^+^^-^"^lyiUf£V',) ) = '-^1' ^° 



^IviUCeVa) ^ 'S'ann-nc {^+^-^\viU{sV2), 

Conversely, if vr does not connect Vi U {eV2) to (— Vi) U (— eV2) and 
'^\v,u(eVo) ^ 'S'ann-nc ( 7+7- ^ I y , jf^Vo J ' Calculate that x(7,7r) = 2. □ 



IViUCeVa) ^ ^ann-nc r+ f- ly^uCeya) 

5.2 Limit distributions 

Asymptotically, the moments of the real ensembles are equal to those of 
their complex analogues. This is not surprising, since intuitively, highest 
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order terms correspond to spheres, which must be orientable, and thus must 
have untwisted edge-identifications, and hence these terms appear in the 
complex expansion as weh. The fohowing lemma gives us an expression for 
the asymptotic value of the moments of any ensemble satisfying ^ (not 
requiring the multiplicativity of /): 

Lemma 5.10. Let {X\}^^^j^ be an ensemble of random matrices satisfying 
with subset of the premaps PMc and function fc- Then for 

7 = (1, . . . ,n) G S'n 

and e : [n] ^ {1, —1} we have: 

Jim E(tr (4f»...4f»)) 

S,n+TzZ'^5^&PMc{±[n]) 

Proof. We know that 

E(tr(x(f»-..4"j")))) = Yl N^^^^'^P'^^-'f,{p). 

pePAM±[n]) 



From Lemma 5.2 x (T) ^epde) < 2, so the exponent on is less than or equal 



to and the limit exists as — )• oo. By Lemma 5.8, the terms which do not 
vanish as — )• oo are those where 5sp5s does not connect [n] and — [n] and 
(5ep(5e|[„] G Iodise- nc (l) ■ Then depSe = tt+ttZ^, and the result follows. □ 

Mixed moments in independent matrices from the various ensembles may 



be expressed similarly using Lemma 4.2 ; however we will not need this result. 



By similar arguments, fluctuations of any matrix ensemble satisfying ([2]), 
including the multiplicativity of /, may be expressed similarly: 

Lemma 5.11. Let {^aIasa '^^ ensemble of random matrices satisfying 
with subset of the premaps PMc and function fc- Let 

7 := (1, . . . , m) (m + 1, . . . , m + n) , 

let 

7op :=(!,..., m) (-m - n, . . . , -m - 1) , 
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and let e : [m + n] — {1, — 1}. Then 



lim A;2 ( TV ( x['^'^^ ■ ■ ■ X^^'^ A , Tr ( • • • xP"'^''^^ 



E 



lim /c (JsTT+vr.^iSe) 

N-^oo 



+ 



E 



TTG S'ann — nc {Top ) 

7r+ vrl (5e e P Afc (± [m+n] ) 



lim /c ((5£7r+7r_^5e) . 



Proof. By Lemma 4.4, we have that 



(e(l)) . . . rp^ / ^(£(m+l)) _ _ _ j^{e{m+n)) 

1 AjTi / ' V + 1 Am+n 



E 



^x{7,5.p4)-4 (p) 



pePAfc(±[m+n]) 
pV { ± [m] ,± [r?i+ 1 ,r?i+n] } = 1 ± [m+„] 

where the second condition under the summation sign means that p must 



connect it [m] and it [m + 1, m + n]. By Lemma 5.2, x (7i ^eP^e) ^ 2, so the 



limit exists as — )• oo. By Lemma 5.9, the terms which do not vanish are 



those where deP^e does not connect ViU{£pV2) and (— Vi)U(— epV2) for some 

choice of sign Sp and 5ep6e\y^^J(^^^v2) either S'ann-nc (t) or 'S'ann-nc (Top)- 

The result follows. □ 

It follows that the three matrix ensembles, or any matrix satisfying ([2]), 
including the multiplicativity of /, has a second-order limit distribution: 

Lemma 5.12. Let {^aIasa be an ensemble of random matrices satisfying 
Q with function fc and subset of the premaps PMc- Then this ensemble 
has second-order limit distribution [A, ipi, ip2) where A is the algebra of non- 
commutative polynomials on indeterminates xx, X £ it A, and ipi and ip2 are 
defined by extending the following expressions by linearity: 



^1 1 4r 



V lim fc{5eTT+TTj6e) 



7I-eS'disc-nc(7) 

6e:7r+7rZ^5.GPAfc(±N) 
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where 7 = (1, . . . , n) and 



E 



7re5 

ann~nc (7) 
<5e 7r+ 7r_ (5e e P Afc (± [m+n] ) 



+ 



E 



where 
and 



— nc (Top) 
(5e 7r+ TT I ^ (5£ e P Afc (± [m+7l] ) 



7 = (1, . . . , m) (m + 1, . . . , m + n) 
7op (1, . . . , m) (-m - 1, ... , -m - n) . 



lim /c {6eTT+7:_^Se) 



Proof. Extending the results of Lemmas 5.10 and 5.11 by linearity, we can 
see that (A, ipi, ip2) meets the conditions on the first two cumulants of traces. 
For r > 3, the cumulant 



kr Tr PI X,,, , . . . , X^,^^^ ,...TT{pr (X,,,,, , . . . , 



is the sum of terms of order x (Tj vr) — r in A^, for various permutations 7 and 



TT. By Lemma 5.2, 7 and vr must satisfy x (7; '''") ^ 2. Thus appears with 
a negative exponent on all terms, which therefore vanish as — >• 00. □ 

Remark 5.13. The moments and fluctuations of these ensembles can be cal- 
culated more explicitly using combinatorial expressions for the number of 
noncrossing diagrams on one or two cycles. In the real Ginibre and GOE 
case, moments and fluctuations are calculated by counting appropriate di- 
agrams. In the Wishart case, terms corresponding to diagrams must be 
weighted by the trace of the matrices along the appropriate permuta- 
tion; however, if we take = Im for all k, the calculation reduces to 
counting appropriate diagrams. 

In each case, the value of the moments is equal to that in the complex 
case, since twisted identifications, or terms in which 5sTr5s connects positive 
and negative numbers, resulting in lower order terms. However, noncrossing 
annular diagrams on both relative orientations of two circles contribute to 
the fiuctuations, and the contribution of each may be different. In the case 
of the GOE and Wishart matrices, the same number of diagrams are possible 
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on both relative orientations, but in the Wishart case transposes of the Dk 
matrices may appear, so the values may be different. 

In the real Ginibre case, the moments are given by counting noncrossing 
*-pairings, that is, those which pair untransposed terms with transposed 
terms. Fluctuations are given by counting annular noncrossing *-pairings 
on the expressions in the two traces as well as those where we have reversed 
the order and the presence or absence of transposes in one of the expressions. 
The number of such *-pairings depends on where the transposes appear. For 
more on *-pairings, see [To] . 

Asymptotically, the GOE has nth moment given by the number of non- 
crossing pairings on n points, that is, for n odd and 0^/2 foi' ^ even, where 
Cfc •= h>fiCk) Catalan number. Fluctuations are given by non- 

crossing annular pairings, which are counted in |21j. Counting both relative 
orientations, the asymptotic value of a fluctuation is given by 



TV 



lim E(Tr(rP),Tr(r'')) 



p+g|,(E_l)!£!(-|_l)| 

4 p! g! 
p+q (E^!)^ (^0^ 
0, 



p, q even 
p, q odd 
otherwise 



The moments of the Wishart ensemble with = Im for all k are given 
by counting noncrossing permutations. These are in one-to-one correspon- 
dence with noncrossing pairings with twice as many points on each circle, 
so the nth moment is C„, and the fluctuations are given by: 



lim K(Tt(WP) ,Tr{W'^)) 



{2p)l (2g)! 



N^CX3 



p + qp\ (p - 1)! q\ {q - 1)1 



6 Preeness 



6.1 First-order freeness 

The real matrix ensembles considered in this paper satisfy the same defini- 
tion of asymptotic freeness as do their complex analogues: 

Lemma 6.1. Let [C] be a set of colours, and associate to each colour c £ C 

a matrix ensemble \Xc^^\ satisfying (2), independent from the other 

I J AeAc l_J 
ensembles. Then these ensembles are asymptotically free. 

Proof. Let : [p] — )• [C] be an alternating word in the colours, and let be 
an element of the algebra generated by the ensemble associated with colour 
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w (k). Then for each /c G [p] we may write 



niH hrik 

i=niA hrtfc_i+l 



for some positive integers ni, . . . , Hp, some A„^-| ^nk-i+ii • • • > '^niH hrife ^ 

A^(j!j), and some function e : [n] — )■ {1, —1} where n := ni + • • • + n^. Let 

4 = [ni H h nfc_i + 1, ni H h n^] and for K C [p], let Ik = \JkeK ^k- 

Let 

7 = 



By Lemma 4.6, we see that 



E tr ii • • • ip = Arx(7,<5.-5.)-2^^ . . . 



TTcePAfJi/,,-! 



(where the last condition under the summation sign means that vr connects 



each block itlfc to at least one other). By Lemma 5.2, we see that this limit 



exists, and by Lemma 5.8, we see that the term associated to vr vanishes 
when N ^ CO unless 5e7r5e/2 G S'disc-nc (t)- 

Consider a vr satisfying the conditions under the summation sign such 
that Ss-kSs G 5'disc-nc (7)- Then (5e7r5e/2 must have a cycle connecting two 
distinct intervals, that is, containing an a and a c such that a G and c £ Ii 
for some k and I, k < I. We know that I ^ k + 1, since w is alternating. The 
permutation (5e7r(5e/2 must have a cycle containing a 6 G Ik+i and a d £ Im 
for some m 7^ + 1. Since <J£'?i"<J£|{a 6 c d} ~ (^)C) (6, d), we must not have 
^\{abcd} ~ ("^5 ^) c, d), SO /c + 1 < 171 < I. By induction, given any connected 
intervals, we can find a pair of connected intervals whose indices are closer 
together, and we derive a contradiction. Thus there are no such vr, and the 
expression vanishes as — )• c«. □ 

The noncrossing conditions satisfied by a term of highest order in an 



expression such as in Lemma 4.6 allow us to further characterize the terms 
contributing to the asymptotic value of the fluctuations, giving us the expres- 
sion in the following theorem. This theorem can be applied to independent 
algebras of any matrices satisfying ([2]), including real Ginibre, GOE, and 
Wishart matrices. 

Theorem 6.2. Let v : [p] ^ [C] and w : [q] ^ [C] be cyclically alternating 
words in a set of colours [C]. To each colour c, associate an independent 
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e 



ensemble of random matrices \Xc'^^\ satisfying (2) with subsets of th 

premaps PMc and function fc- For integers mi, • • • , mp and rii, • • • , Ug, let 
m := mi + • • • + mp and n := ni + ■ ■ ■ + Uq. Let e : [m + n] — t- {1, — 1}. 
For each A: G [p] and i E [mi + • • • + mk-i + 1, mi + . . . + m^] let Aj G ^v(k) 
and for k G [q\ and i £ [m + ni + ■ ■ ■ + rij^^i + 1, m + ni + • • • + n^] let Aj G 

miH hmfe 

= n 

i=miH hmfc_i+l 

and 

m+niH hrifc 

n 

i=m+niH |-nj,_i+l 

r/ien if p ^ q 



lirn^ k2 [Tr \^Ai---Apj,Tr[Bi--- Bg 
and if p = q > 2, 

lirn^ k2 (Tt (ii • • • ip) , Tr • • • Bp 

p-i p 

VTT lim (E(tr(^,Bfc„i))-E(tr(A,))IE(tr(i?,,_,))) 



■ 7V-5-00 
fc=Oj=l 



+ E n - E (tr (^.)) K (tr {Bl^,))) . 

^ — ' _/v—>oo 
fc=Oi=l 

Proof. For integer /c (taken modulo p), let 

Ik ■■=[mi-\ h mfc-i + 1, mi H h m^] 

and for k (taken modulo q), let 

Jk ■■= [m + ni -\ h nfc-i + 1, m + ni H h nfe] . 

For -ftT C [p], let Ik = UkeK^k, and for K C [q], let Jx = \JkeK ^K- 
Let 

7 = (1, . . . , m) (m + 1, m + n) . 
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By Lemma 4.6, we know that 



k2 [Tt (^^1 • • • Apj , Tr (^Bi • • • 

^cePA/e(/„_i(^)UJ„_i(^)) 

7rV{±H.±[™ + l.™+"]} = l±[m + n] 

±/fe,±jfc^7rv{±/af^,u{±jaf^i 

(the third hne under the summation sign says that vr connects the blocks 
± [m] and it [m + 1 , m + n] , and the fourth says that vr connects each block 



zblfc and ibJfc to at least one other). By Lemma 5.9 we know that the terms 
that survive in the limit N ^ oo are those such that (5e7r(5e doesn't connect 
[m] U (Ett [m + 1, . . . , m + n]) to — [m] U (— [m + 1, m + n]) for some choice 

of sign e^, and dsirds G 5'ann-nc ( 7+7^^ I 



[m]U(£-n-[m+l,m,+n])^ 

We now wish to show that for such a vr, JeVr^g must connect intervals 
in "spokes" as shown in Figure [T| where each spoke can be expanded to 
a noncrossing diagram which connects the those two intervals. Let ■ = 

^+^- l[m]U(£^[m+l,m+7i])' 

First we show that 5sTt5s cannot have an orbit connecting two distinct 
Ik or two distinct e^Jk- By the arguments in Lemma |6.1[ we can find an 
a £ Ik, b £ Ik+i, c £ Ii and d such that (^eTT^elja & c d} ~ (o, c)(6, d), and 
here we must have 77r|{a bed} ~ (^i ^' ^) (^)- Reversing the roles of a and c, 
we may then find an x S J^+i sharing an orbit of 5£TT~^5e with a y from a 
different interval, and we conclude similarly that 77r|{a,e,x,y} 
Since 7n\{a,b,c,x} = ia,b,c,x), Ss7r6s\{a,b,c,x} = (a'C)(6)(x) (first annular- 
noncrossing condition), so dgiTdel^^ bcdxy} ~ ^) f)- other 

hand, Ax,y||^ bed} ~ ^' '^)' violating the third annular- noncrossing con- 
dition. So there must be no such orbit. 

We now show that SeTtSe may not connect an interval in one orbit of 
77r to two distinct intervals in the another. If so, let c and y belong to the 
same interval (in the orbit of which we will call 7ext) and share cycles of 
SsTt6s with a and x respectively, where a and x belong to distinct intervals 
in the other cycle of which we will call 7int (thus these cycles of 6s'n'Ss 
must be distinct), and such that if we apply 77r repeatedly to y we will find 
c before we find an element of another interval. If we apply 7,^ repeatedly 
to a, we must find various b belonging to another colour before we find x. 
At least one of these b must be connected by a cycle of SsTtSs to another 
interval, which by the above arguments must be in 7ext, and because it is 
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of a different colour from a, c, x and y, in an interval distinct from the one 
containing c and y. We have thus that SeirSel^^ bcdxy} ~ i'^^^) (^j ^) u)- 
However, by our arguments, ^x,y\[ab cd} ~ i^^^i'^jd), violating the third 
annular-noncrossing condition. 

We find thus that each interval in an orbit must be connected by S^irS^ 
to exactly one other, which must be in the other orbit of 77r. We conclude 
that p = q for any nonvanishing covariance. 

We now show that the diagram on the intervals will be a spoke dia- 
gram: if X G /fc and y G Et^Ji share an orbit of SsTtSs, then I^^i must be 
connected to e^J/_e^. Assume not: then SsTtSs must connect Ik+i to an- 
other e-j^Ji' and e^J^_£^ to another I/./. Let a G Ik+i and c G EnJi' share 
an orbit of SsTtSs, and let d G e-yrJl-e-n ^ ^ ^k' share another. Then 
<^^^^£l{a,fe,c,d,x,?;} = (O'C) (a;,y). However, A^,j,|^^ ^^^^^^^ = (a,6,c,d), vio- 
lating the third annular-noncrossing condition. 

We now show that each spoke consists of a noncrossing diagram on and 
e-i^Ji. The permutation -^x,2/|/^u(£^jj) is identical for any choice of x G 7ext 
and y G 7int not in /fcU(e,r^O' since we are considering^ = q >2, such an 
X and y will exist. We will thus refer to this permutation without specifying 
X and y. The second and third annular-noncrossing conditions become the 
disc-standard and disc-noncrossing conditions relative to Xrvlr t\- 

Conversely, we show that any premap vr such that SsTrSel^^ has such a 
spoke arrangement for some choice of sign (that is, one which connects li 
to only s-n-Jk-cni ■"*ce versa for some k and such that Se7TSs\j.u^^^ji__^ .) G 

'S'disc-nc (A^,2/l/iU(£^Jfc_e„i)) ^ ^ ^^ext and y G 7int with x,y ^ li U 

{e^^Jk-e^i)) must be in S'ann-nc (77r)- 

Any restriction of \x,y to Ik or e-j^Ji is equal to 77r restricted to the same 
domain. Any ^eTr^e satisfying the disc-standard condition on X^^y will then 
satisfy the first annular-standard condition. 

If a, 6, c and d share a cycle of 5e7r(5e, with a and 6 in one cycle of 7,r 
and c and d in the other, then a and 6 must be encountered in some order 
in Xx^y before any elements of the cycle of containing c and d, so either 
(a, 6, c) or (a, d, 6) is disc nonstandard on Xx,y Thus (JeTr^e must satisfy the 
second annular-standard condition. 

The elements of any spoke in one of the cycles of 77r comprise an interval 
of that cycle, so if 77r|{a,6,c,d} = ia,b,c,d) and SeTrS^l^^^^^^^d} = (a,c)(6,d), 
then one of b and d (and hence the other) must be in the same spoke as a 
and c. The elements of this spoke must be disc-crossing relative to Xx,y, so 
S^ndg must satisfy the first annular-noncrossing condition. 

If x' and y' sharing a cycle of devr^g are on different cycles of 77r, they 
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divide the other elements of their spoke into two intervals in the cycle 
'^^.yl/feU(e Ji)' Since JeVrJe is disc-noncrossing relative to 

Xx,y, any cycle of iJeTnJg other than the one containing x' and y' must be 
contained in one of the Kj. We note that each Kj intersects each cycle of 
7^ in an interval not containing x, y, x' or y' , so both Xx^y and \x',y' induce 
a permutation on Ki in which elements are mapped under 7,^ within the 
interval belonging to a cycle of 7,^ then the last element of an interval is 
mapped to the first of the other. Thus Xx',y'\j^. = ^x,y\x.- Since 6^716^ 
must then be disc standard relative to Xx,y for any x and y, even those in 
the same spoke, it must satisfy the second annular-noncrossing condition. 

We note that in the cycle of Xx,y, where x £ li and y G s-wJk-e-nii we 
see elements from the intervals in the cyclic order /j, /j+i, . . . , h-i, /«, Jk-e^^i, 
Jfc_e^(j_i), . . . , Jk-e^{i+i), Jk-e-ni (recalling that subscripts are taken modulo 
their appropriate range) . If there is an instance of the third annular-crossing 
condition (variables as given in the definition) where a and c belong to a 
different spoke from x and y, then h and d must belong to the same spoke 
as a and c: if a and c belong to the same interval of their spoke in Ax,j/, then 
one of 6 or d must appear in that interval; and if a and c belong to the two 
different intervals, then no two elements of any other spoke may appear in 
the correct order. In either case, they are disc-crossing relative to Xx^y, so 
this configuration cannot occur. If a and c belong to the same spoke as x 
and y, then they must be contained in one of the Kj^ as above. We note 
that the Kj coincide with the intervals of this spoke in Xx^y, so b and d must 
also be contained within this interval. However, if 6^715^ is crossing relative 
to Ax,y, it is crossing relative to Xx'^y' for x' and y' in another spoke, so this 
configuration also cannot occur. 

We now rearrange our expression for the asymptotic covariance in terms 
of the spokes. Since a vr contributing to the asymptotic value of the co- 
variance must connect [m] to exactly one of [m + l,n\ and — [m -|- l,n], tt 
uniquely determines the sign e-,^ and the value k such that each /j is con- 
nected to ej^Jk-ET^i for all i G [p], and thus the asymptotic covariance can be 
expressed as a sum over all spoke diagrams for each choice of sign and k. 

For a given eq and A;, the contribution is a sum over terms corresponding 
to a choice of connected spoke on each Ij and e^Jk-eoi (where the contribu- 
tions are multiplied) , and as such may be factored into sums of all connected 
spokes on Ii and e^Jk-eoi fo^ all z G [p]. 

A disc-noncrossing permutation on 1^ and EqJi which does not connect 
the two intervals induces disc-noncrossing permutations on each of Xx,y\j^ 



and Xx,y\i,^j^- In terms of the expression given in Lemma 



5.10 



the sum over 
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those which do connect the two intervals is then 



hm E tr L4feR 



(eo) 



Hm E(tr(^A:)) hm E tr LB} 



where Lemma 5.10 guarantees the existence of this limit if v (k) = w (l) and 
it is equal to zero for all N if v (k) ^ w (l). The result follows. □ 



7 The main definitions 



We take the property satisfied in Theorem 6.2 as our definition of asymptotic 
real second-order freeness: 

Definition 7.1. Let \ Xc^^ \ be an ensemble of random N x N matri- 
l J AeAc 

ces for each colour c G [C]. We say that the ensembles are asymptotically 
real second-order free if they are asymptotically free, have a second-order 
limit distribution, and if, for any v : [p] ^ [C] and w : [q] ^ [C] cycli- 
cally alternating words (or words of length 1) in the set of colours [C] and 
Ai, . . . , Ap and Bi, . . . , Bg random matrices with A^ in the algebra gener- 
ated by the ensemble associated with v (k) and B^ in the algebra generated 
by the ensemble associated with w (k), we have for p ^ q, 

llmj, {Tr[A,...Ap),Tr[B,...Bg))=0 

and for p = q > 2, 

limj, {Tt[A,...Ap),Tr{B,...Bp)) 

= E n J™ ^ {A,Bk-i)) - E (tr {Ai)) E (tr {Bk-^)) 

I — n ^'—1 \ / 



fc=0 j=l 

+ E n f J™ ^ {A,Bl^,)) - E (tr {A,)) E (tr (5,^+,)) 

fc=Oi=l ^ 

Noting that 

^ _ Iff w 4 . ( R(±l)^ ^ _ r j . r(±i) 



E [ii [AkBl^^'"^"^ - E (tr {Ak)) E (^tr j j = ^ [^kB\ 

the above condition is equivalent to the following condition on the algebra 
generated by the second-order limit distributions of the matrices, which we 
take as our definition of second-order freeness: 
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Definition 7.2. Let Ai,.. . , Ac be subalgebras of A, {A, ipi, ip-i) a second- 
order noncommutative probability space equipped with an involution a i— t- a* 
reversing the order of multiplication. Then ^i, . . . , Ac are real second-order 
free if they are free and if, for every ai, . . . ,ap,bi, . . . ,bq G A such that 
«fe € and bk G A^(fc) for v : \p] ^ [C] and u; : [g] ^ [C] cyclically 

alternating words (or words of length 1) in [C], for p q, 

cp2 (ai--- hp, bf-bq^ =0 

and for p = g > 2 

1 p p^i p 

</^2 (ai--- hp, bi---bp^ ^'^Yi^^ (aih-ij + ^Y['Pi ■ 

k=0 i=l k=0 i=l 
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